Nonanalytic Magnetic Response of Fermi- and non-Fermi Liquids 
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We study the nonanalytic behavior of the static spin susceptibility of 2D fermions as a function of 
temperature and magnetic field. For a generic Fermi liquid, X f(T, H) = const + ci max{T, /j,b \H\}, 
where c\ is shown to be expressed via complicated combinations of the Landau parameters, rather 
than via the backscattering amplitude, contrary to the case of the specific heat. Near a ferromagnetic 
quantum critical point, the field dependence acquires a universal form xJ 1 (^) = const — C2\H\ 3 ^ 2 , 
with C2 > 0. This behavior implies a first-order transition into a ferromagnetic state. We establish 
a criterion for such a transition to win over the transition into an incommensurate phase. 

PACS numbers: 71.10. Ay, 71.10 Pm 
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The nonanalytic behavior of thermodynamic quanti- 
ties of a Fermi Liquid (FL) has attracted a substantial 
interest over the last few years. The Landau Fermi- 
liquid theory states that the specific heat coefficient 
7 (T) = C(T)/T and uniform spin susceptibility X s(T, H) 
of an interacting fermionic system approach finite values 
at T, H = 0, as in a Fermi gas. However, the temperature 
and magnetic field dependences of j(T, H) and Xs(T, H) 
turn out to be nonanalytic. In two dimensions (2D), 
bothj and Xs are linear rather then quadratic in T and 
\H\ [lj. In addition, the nonuniform spin susceptibility, 
Xs (q) , depends on the momentum as \q\ for q — ► 

Nonanalytic terms in 7 and Xs arise due to a long-range 
interaction between quasiparticles mediated by virtual 
particle-hole pairs. A long-range interaction is present in 
a Fermi liquid due to Landau damping at small momen- 
tum transfers and dynamic Kohn anomaly at momen- 
tum transfers near 2kp (the corresponding effective in- 
teractions in 2D are |0|/r and |f2| cos(2kpr)/r 1 / 2 , respec- 
tively). The range of this interaction is determined by the 
characteristic size of the pair, L p h, which is large at small 
energy scales. To second order in the bare interaction, 
the contribution to the free energy density from the inter- 
action of two quasiparticles via a single particle-hole pair 
can be estimated in 2D as SF ~ u 2 T/L^ h , where u is the 
dimensionless coupling constant. As L p h ~ vf/T by the 
uncertainty principle, SF oc T 3 and 7 (T) oc T . Likewise, 
at T = but in a finite field a characteristic energy scale 
is the Zeeman splitting /ib\H\ and L p h ~ Vp / jj,s\H\. 
Hence 5F oc \H\ 3 and X s (H) oc \H\. 

A second-order calculation indeed shows 0, 0, Bl that 
7 and Xs do depend linearly on T and \H\. Moreover, 
the prefactors are expressed only via two Fourier compo- 
nents of the bare interaction [f (0) and U(2kp)] which, 
to this order, determine the charge and spin components 
of the backscattering amplitude T CiS (9 = w), where 9 is 
the angle between the incoming momenta. Specifically, 
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where 5j(T, H) = j(T,H) - 7 (0,0), S X s(T,H) = 
Xs(T,H) - Xs (0,0), 6 Xs (q) = Xs(q) - Xs(Q), r c (n) = 
(m/2n) [U(0) - U{2k F )/2], T s (tt) = -(m/47r)C/(2fc F ), 
7(0) = to7t/3, Xs(0) = n 2 B m/ir,E F = mv%/2, \i B 
is the Bohr magneton, and the limiting forms of the 
scaling functions are / 7 (0) = f x (0) = 1 and / 7 (x ^> 
1) = 1/3, f x (x ^> 1) = 2x. (Regular renormaliza- 
tions of the effective mass and g— factor are absorbed 
into 7 (0) and Xs(0)). The second-order susceptibility in- 
creases with H and q, indicating a tendency either to a 
metamagnetic/first-order ferromagnetic transition or to 
a transition into a spiral state. These tendencies signal 
a possible breakdown of the Hertz-Millis scenario of the 
ferromagnetic QCP 

Experimentally, a linear T dependence of the specific 
heat coefficient has been observed in thin films of 3 He 0- 
A linear increase of Xs with magnetization (and thus 
H) has been observed in a 2D GaAs heterostructure 0. 
Since none of these experiments correspond to the weak- 
coupling limit, there is obviously a need for a nonpertur- 
bative treatment of nonanalytic terms. 

It has recently been shown [J, |9| that the second-order 
result for j(T, 0) in Eq.Q becomes exact once the weak- 
coupling backscattering amplitudes r c , s (7r) are replaced 
by the exact ones. This implies that the 0(T) term in 7 is 
determined exclusively by ID scattering events embedded 
into the 2D space. In these events, fermions with almost 
opposite momenta experience either a small or almost 
2kp momentum transfer. It has been conjectured in Ref. 
|j| that the nonanalytic part of the spin susceptibility can 
be generalized in the same way, i.e., by replacing weak- 
coupling T s in Eq. Q by its exact value. The same 
result was obtained within the supersymmetric theory of 
the Fermi liquid [lfj , and in the analysis of the scattering 
amplitude in the Cooper chanel, T^(0) 11]. Since an 
extension of the second-order result for 5 Xs hints at far- 
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reaching consequences for the ferromagnetic QCP, it is 
important to establish a general form of Xs {T, H) . 

In this communication, we present a general analysis 
of the nonanalytic behavior of the spin susceptibility in 
2D. We show that, in contrast to the specific heat case, 
higher orders in the interaction are not absorbed into the 
renormalization of r s (7r) (equal to Ff (it)), but give rise 
to extra 0(T, H) terms, whose prefactors are given by in- 
finite series of the scattering amplitudes F s (9) and F^(6) 
at all angles. That higher-order terms in F~(9) are rele- 
vant was noticed in Ref.|ll|. We show that higher-order 
terms in F s (0) are also present. These terms are more 
important than higher powers of F^(9) as the latter are 
logarithmically reduced at small T and H. When the 
interaction is neither weak nor peaked in some particu- 
lar channel, the total prefactor of the 0(T, H) term may, 
generally speaking, be of either sign. However, the uni- 
versality is restored near a ferromagnetic QCP, where the 
n = partial component of F s {6) diverges. We show that 
near the QCP the inverse susceptibility X^ 1 ^ = 0,-ff) 
behaves as £~ 2 — A|i7| 3 / 2 , where £ is the correlation 
length. This dependence is dual to the £ -2 — B \q^ 2 be- 
havior of the nonuniform susceptibility Signs of A 
and B are positive, so that the nonanalytic terms destroy 
a continuous transition towards a uniform ferromagnetic 
state, and, depending on parameters, the system under- 
goes cither a second-order transition into a spiral phase 
or a first-order transition into a ferromagnetic state. 

The temperature and magnetic-field dependences of 
j(T,H) and Xs(T,H) are most straightforwardly ob- 
tained by evaluating the thermodynamic potential at fi- 
nite T and H, 5 (T, H) , and then differentiating it twice 
with respect to T or H, respectively. To understand the 
difference between the spin susceptibility and the spe- 
cific heat, consider for a moment the case of a contact 
interaction: U(q) = U. To second order in U, the ther- 
modynamic potential 5 is expressed via the convolutions 
of the polarization bubbles II (fi m , q, H) (with opposite 
spin projections for H ^ 0). The polarization bubble has 
a conventional form 
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For large momenta (vf<7 3> |fi m | ~ ^b\H\), the dy- 
namic part Ildyn behaves as |r2 m |/g. Consequently, the 
momentum integral J d 2 <7ll 2 yn diverges logarithmically 
and is cut at q = max{\Cl m \, fi B \H\}. Because of the 
logarithm, the subsequent summation over frequencies 
yields a universal term 5(T, if) cx max{T 3 , (/j, B \H\) 3 }. 
More precisely, 3{T,H) cx T 3 p{^ B H/T), where p(x < 
1) = a + bx 2 + . . . and p(x ^> 1) ~ |x| 3 . Accordingly, 
<5 7 (T)/ 7 (0) ~ 8xs{T)/xs(0) cx T, and 5 Xs (H) cx \H\. 



To second order, Sj(T) and 5xs(T,H) behave simi- 
larly. The difference between the two quantities shows up 
at higher orders in U. In the rest of the paper, we con- 
sider only scattering in the particle-hole channel. As we 
have already mentioned, there are higher-order contribu- 
tions from the Cooper channel, but they are logarithmi- 
cally small in a generic Fermi liquid, and nonsingular near 
a ferromagnetic instability. A particle-hole contribution 
of order n contains integrals of II™ = Ilg + d-ill^ 1 + 
where c n are the constants. The nonanalytic part of 
7(T) is solely related to the logarithmic divergence of the 
momentum integral J d 2 qU 2 yn cx log|fi m |, because only 
the log singularity ensures the nonanalytic result of the 
subsequent Matsubara sum: T fi^ log |fi m | =const- 
0(T 3 ). The momentum integrals of II*j yn with k > 2 
are not logarithmically divergent, and the subsequent fre- 
quency summation gives rise only to regular, T 2 ,T 4 ... 
powers of T in j(T). As a result, higher-order diagrams 
for 7(T) only renormalize bare interaction U into a full 
backscattering amplitude. For Xs(H), the situation is dif- 
ferent - higher powers of Ildyn do contribute additional 
|iJ| 3 terms to S, and hence additional \H\ terms to the 
susceptibility. Indeed, evaluating 3(0, H) to third order 
in U and retaining only the contribution <5S( 3 ' 3 ) (0, H) 
from n| yn , we obtain 
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where u = mU/(2ir). The momentum and frequency 
integrals in 5S^ 3,3 ^ come from the region |f2 m | ~ vpq ~ 
fi B \H\. This implies that the U 3 \H\ 3 term in J3J appears 
by purely dimensional reasons, and does not require the 
g-integral in 2 to be logarithmically divergent. 

We see that the \H\ terms in Sxs(H) coming from two 
and three or more dynamic bubbles correspond to phys- 
ically distinct processes. The distinction becomes im- 
portant for a generic Fermi liquid. The contribution to 
Sxs{H) from two dynamic bubbles, which starts at or- 
der U 2 , is generalized beyond the weak-coupling limit 
by replacing the bare interaction by a fully renormalized 
static vertex. Using the same computational procedure 
as in Ref.0], we find that, similarly to the specific heat, 
the exact result for this contribution is expressed in terms 
of the spin part of the backscattering amplitude F s (n) : 
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The same result was obtained in Refs. The log- 

arithmic divergence of the momentum integral of II^ yn 
is the crucial element in the derivation of Eq.(0}, as 
the higher-order corrections can be absorbed into static 
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r s (7r) only if typical v F q are much larger than typical 
|fi OT |, given by hb\H\. 

On the other hand, contributions to 8xs(H) from three 
and more dynamic bubbles come from v F q ~ |^m|> an d 
are expressed via the convolutions of the partial harmon- 
ics of the scattering amplitudes, which do not reduce to 
higher powers of the backscattering amplitude. As an il- 
lustration, we consider a generalization of the third-order 
contribution 5!E/ 3 ' 3 ', assuming that the spin component 
of the scattering amplitude has only two partial compo- 
nents: n = 0,1, i.e., T s (8) = r s , + cos#r Si i. Replacing 
each interaction vertex by T s (9), we obtain 
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x [r 3 i0 - aiiW - a 2 r s , r 2 x + a 3 Tl J . (5) 

where a\ = 3(21n2-l), a 2 = 3(31n2 - 2), a 3 = 
(5/2 — 3 In 2). This expression obviously does not reduce 
to the cube of the backscattering amplitude, which in this 
approximation would be equal to r s (-7r) = r Si o — r s ,i- 
Higher-order contributions are given by progressively 
more complicated combinations of T Si o and r S) i. 

The total result for 8x s is a sum of backscattering 
and all-angle scattering contributions. Since r s (-7r) is 
equivalent to the spin component of the particle-particle 
scattering amplitude T%(tv) (Refs. 0, the repul- 

sive interaction in the Cooper channel leads to a log- 
arithmic reduction of T s (tt) [H El El For T = 

0, if -> 0, TJn) cx l/log|if| and therefore Sxf S oc 
| if | / log 2 |if| [l4(- On the other hand, contributions to 
5xs from three and more dynamic bubbles contain an- 
gular averages of T s (9) , which are not affected by the 
Cooper singularity. Therefore, the |if | terms from these 
contributions do not acquire additional logarithms and 
win over the backscattering contribution for T, if — > 
(and also over higher-order terms in T^(9), which vanish 
logarithmically at T, if -> 0). Note that for |r s ,i/r s>0 | < 

1, the sign of xf ny in Eq. JHJ) is determined just by the sign 
of i\o : for negative r Si o (corresponding to enhanced fer- 
romagnetic fluctuations), Sx s increases with if, whereas 
for positive r Si o, $Xs decreases with if. 

Next, we discuss the behavior of the spin susceptibil- 
ity in the vicinity of a ferromagnetic QCP, where r Sj o 
diverges, while other components of T s remain hnite. At 
any hnite distance from the QCP, the backscattering am- 
plitude still vanishes as l/logmax{^s|if |,T}. However, 
at large r s o, this behavior is confined to an exponen- 
tially small range of if and T, which we will not con- 
sider below. Outside this range, the backscattering am- 
plitude diverges as r s ,o, and the backscattering contri- 
bution Sxf S (T,H) diverges as T 2 . All- angle contribu- 
tions, however, diverge even stronger, and one needs to 
sum up a full series of diagrams to obtain the behavior 
of 6Xs(T, if) near QCP. To do this, we assume, as it was 
done in Refs. [3. Il5j. that the Eliashberg approximation 



is valid near QCP because overdamped spin fluctuations 
are slow compared to fermions. In the Eliashberg theory, 
the field-dependent part of the thermodynamic potential 
is H = (l/27r)T£ n fdgg log Xo'Oz.fi.-fO , where 



Xo(q,ft m ,H) = 



tt S + a 2 q 2 + (27r/m)IIdy n (fi m , q) ' 



(6) 



is the dynamic spin susceptibility without nonanalytic 
corrections, 8 = |r s .o| _1 , and a is the radius of the ex- 
change interaction, required to be large (akp 3> 1) for 
the Eliashberg theory to work 16] . nd yn differs from 
Eq.(J2J in that (i) it is built on full Green's functions 
(containing self-energies) and (ii) /is in the denominator 
is replaced by fi* B — fi B /5 (cf. Ref. 01). To illustrate 
once again the difference between the specific heat and 
spin susceptibility, we set £ = oo in the denominator of 
Eq.@) and neglect the self-energy renormalization for a 
moment. Evaluating the derivatives of S(T, if) with re- 
spect to T and if, we find that the prefactor of the T 
term in the specific heat coefficient diverges 0, whereas 
the prefactor of the |ff| term in the spin susceptibility 
remains finite: 5xs(H)/(^* B f = (2/nv 2 F )\H\. This indi- 
cates dramatic cancellations between diverging terms in 
the perturbation theory for Sxs(H) |is|. 

A complete result for the susceptibility near QCP is 
obtained by including the self-energies when evaluating 
riny n in Eq.© (vertex corrections are small, see Ref. 
[l5j). The self-energy near the QCP interpolates between 
£ = Xu> m away from QCP andS = WQ^ 3 |w m | 2 / 3 near 
QCRwhere A = 3/(Ak F aS) and lu q = 'SV3E F / '(A(k F a) 4 ) 
|l9i |2PJ. Using these expressions, we obtain for the in- 
verse susceptibility 



' 1 (H,T = 0) cx 8 



3 v F /a 



where K H (0) = 1 and K H (x > 1) = 1.25y/x, The limit 
of x — * oo describes the situation right at QCP. Here, 
divergent £ cancels out from the answer, and the if de- 
pendence of xj 1 becomes |ff | 3 / 2 . We emphasize that the 
exponent of 3/2 is the consequence of the non- Fermi liq- 
uid behavior, manifested by the divergence of the "effec- 

— 1/3 

tive mass" dY,/du> m cx u> m 

The nonanalytic |if | 3 ^ 2 dependence exists only at T — > 
0. At finite T, the field dependence of the spin suscep- 
tibility is analytic: 8x s oc if 2 . However, the prefactor 
scales as I /(XT) away from QCP, and as T" 1 / 6 at QCP. 
AT if = 0, 8 Xs (T) oc TlogT [Hi!. 

A complementary way to see the nonanalytic depen- 
dence the susceptibility on the magnetic field is to an- 
alyze the thermodynamic potential itself. Viewed as a 
function of the magnetization M — (m/7r)rj/XB, where 
2r\ is the difference of the Fermi energies for spin-up 
and spin-down fermions, the thermodynamic potential 
S(T = 0,?y) contains a nonanalytic |?y| 3 term away from 
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criticality: S nQ (0, rj) = -\r)\ 3 /4Snv F X. Near QCP, A di- 
verges and the |7y| 3 dependence is replaced by |ry| 7 ^ 2 , in 
agreement with the H 3 ! 2 field dependence of the spin 
susceptibility. At finite T, the \r]\ 7 / 2 term evolves into an 
analytic rp 4 one with a singular prefactor T -1 / 6 . 

We now study the consequences of the nonanalytic be- 
havior of Xs(T, H, q). First, we see from Eq.0 that the 
spin susceptibility diverges at some finite value of H, 
which implies that a second-order ferromagnetic QCP is 
preempted by the first-order one. This possibility was 
discussed in detail in Ref.[2^] - our analysis differs from 
this work in that we include the fermionic self-energy and 
nonanalytic T dependence of the susceptibility. Assum- 
ing that the first-order transition occurs near QCP, where 
the nonanalytic term in S(T = 0, 77) is rp 1 1 2 , we have 

^=(T = 0, ,) = |,»-!§^ + 6V, (8) 

with E — 3.82Ep / (kpa) 4 ^ 3 ■ Because of the nonanalytic 
term, S has a minimum at finite 77. The first-order tran- 
sition occurs at 5h = (fcpa/3.32) 8 /(bEp) 6 when 5 = 
at this minimum. By order of magnitude, b ~ 1/Ep. 
The first-order transition occurs in the critical region 
(S H < 1) for k F a < 3.32(bE F ) 3 / 4 . 

The susceptibility at H — but finite q is given by 

X: 1 {q)=5 + a 2 (q 2 -cq 3 / 2 k]l 2 ) 1 (9) 

where c ss 0.25 15]. Xj 1 ^) diverges at q = qo = 0.035fc_F 
for 8 q ~ 0.42 x 10~ 3 (akp) 2 ■ This signals a transition into 
an incommensurate phase. 

Which of the two instabilities occurs first de- 
pends on the nonuniversal parameter p = 5 q /8u = 
(l.35bEp I 'akp) 6 . For p > 1, the first instability is into 
the incommensurate phase; for p < 1, the first-order 
transition occurs first (see Fig. QJ. Although formally 
akF should be large, both situations are actually possi- 
ble, particularly if bE F is a large number. 

At finite T, the transition line has an S— shaped form 
(see Fig ; _|3 because of the negative T dependence of 
X _1 (T) I23. The tricritical point separating the second- 
and first-order transitions, results from the balance be- 
tween the b 2 r\ 4 term in Eq. (|SJ) and the rf'/T 1 / 6 term 
which replaces the \r]\ 7 ^ 2 term at a finite T. 

To summarize, in this paper we considered the tem- 
perature and magnetic field behavior of the spin suscep- 
tibility of a 2D Fermi liquid, both away and near a fer- 
romagnetic QCP. We found that in a Fermi-liquid phase, 
Sxs(T,H) oc max{T, \H\\], but the prefactor is not ex- 
pressed solely in terms of the backscattering amplitude, 
in contrast to the specific heat. At a ferromagnetic QCP, 
the magnetic field dependence of X _1 (T = 0, H) becomes 
H 3 ^ 2 , with the universal, negative prefactor. This behav- 
ior favors a first-order transition to fcrromagnetism and 
competes with q 3 ! 2 behavior of x~ X {T = H = 0,q) which 
favors an incommensurate spin ordering. 




FIG. 1: Schematic phase diagram near a ferromagnetic QCP 
for p < 1 (see text). 8 is a control parameter, e.g., pressure. 
The second-order phase transition becomes first-order below 
a tricritical point because of the nonanalytic r] 7 ^ 2 term in 
Eq. (|HJ . At small T, the transition line has an S'-shaped form 
because of the negative T dependence of X^ 1 (r). 
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